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Abstract. Gravitational lensing will cause a dispersion in the Hubble diagram for high redshift sources. This effect 
will introduce a bias in the cosmological parameter determination using the distance-redshift relation for Type la 
supernovae. In this note we show how one can diagnose and correct for this bias when doing precision cosmology 
with supernovae. 

Key words. Gravitational lensing - cosmological parameters - Methods: statistical 



1. Introduction 

During the last decade, gravitational lensing has become 
one of the most important tools in cosmology. Weak lens- 
ing measurements are used to obtain information on the 
amo unt and distribution of mat t er in the univer s e (see , 
e.g., van Waerbeke et al., 2001; Hennawi et al, 2001). 



Strong lensing has long been recognized as a tool to 



constrain cosmological parameters (see, e.g., Koopmans 



fc Fassnacht, 1999|; Browne, 2000, for measurements of 



the Hubble constant, Hq). Controversial upper limits on 
the amount of dark energy have been derived from the 
statistics of multiply imaged quasars and radio sources 



' Kochaneck , 


1996; 


Falco et al. 


1998) 



set of Schechtcr parameters for E/SO galaxies Chiba & 
|Yoshii (199E) conclude that the same quasar data set 
is best fitted with a universe dominated by dark en- 
ergy, i.e., in agreement with the concordance cosmology, 
{nM,^x) = (0.3,0.7) whereas |Keeton| ( |2002D argue that 
when calibrated to counts of distant galaxies, the method 
loses most of its sensitivity to the amount of dark energy. 

Further, in Goobar et al. (20021) it is shown how one 
could use multiple images of supernovae (SNe) to simul- 
taneously measure Hq, the fractional energy density com- 
ponents and the equation of state parameter of a possi- 
ble dark energy component. The classical SN cosmological 
test however, is to study the distance-redshift relation us- 
ing Type la SNe as standard candles. The effect of gravita- 
tional lensing on these measurements is to cause an addi- 
tional dispersion in the observed magnitudes and thus be 
a source of systematic error in the cosmological parameter 



et al 



1997; Mortsell 



determination (see e.g. Wambsganss et al 

2001b| ). In this note we show how the negligence 



of this effect may lead to a systematic underestimation of 
JIm- We also investigate how this bias can be diagnosed 



and corrected for by using the correct probability distribu- 
tion function (pdf) for the dispersion in SN magnitudes. 

In Sects. ^ and |^ we discuss the method and the nu- 
merical simulations used in the analysis. In Sect. ^, we 
investigate what kind of bias we expect if not taking lens- 
ing effects into consideration, i.e., assuming a Gaussian 
magnitude distribution. In Sect. ^, the dispersion in mag- 
nitudes due to gravitational lensing is parametrized for 
different matter distributions and redshifts and in Sect. |^ 
we correct for lensing effects by performing a maximum 
likelihood (ml) analysis using the correct pdf's. 



2. The distance-redshift cosmological test 

The distances and redshifts of Type la SNe can be used 
to fit, e.g., the mass energy density, ^Im, and the dark en- 
ergy density, fix, in a Friedmann-Lemaitre universe. The 
relation between the parameters can be expressed as 



m{e,M,z) = X + 51ogioK(0,2)] , 
M = 25-FM + 51ogio(c/i?o) 



(1) 
(2) 



where is a nuisance parameter containing the absolute 
magnitude, M, of the SNe and Hq is the Hubble constant. 
Here 6 represents the parameter vector 6 = {Qm, ^x, w). 
Further the luminosity distance d'j^ is given by 



d'L = 

/ = 



(l + z)^J=sin(^/^/) 
(1 + ^)/, 

(l + z)^sinh(Vnr/), 

1 — r^M — , 

dz' 



nk < 
nk = o 
nk>o 



(3) 

(4) 
(5) 
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H\z) = ^(1 + z)3 r!M + nx{i + z)3(-+i) + (1 + nk . 
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The usual way to proceed is to make a x fit of 
{ilM,^x,-M) using Eq. (Q) assuming that the dispersion 
in magnitudes is Gaussian distributed (as the intrinsic dis- 
persion of Type la SNe seems to be). The same procedure 
can be used to fit, e.g., the equation of state parameter. 



w, of a, more general dark energy component (Goliath el 
al 



20(1) 



Gravitational lensing will however induce an asymme- 
try in the magnitude distribution and in order to avoid a 
bias in the cosmological parameter estimation, we need to 
perform a ML-analysis with the correct pdf. 



to the predicted results from the Supernova Factory 
Campaign ( Aldering , Supernova Factory Webpage ) that is 
scheduled to start in 2003. Since lensing effects get larger 
at higher redshifts, the magnitudes of theese SNe are un- 
affected. For example, for a supernova at redshift x = 0.08 
the lensing effects are in general of the order of 0.01 mag- 
nitudes for a dark matter model where the fraction of 
compact objects is assumed to be 20%. 

All simulations have been made assuming a universe 
with = 0.3, Q,x = 0.7, w = —1 and Hq = 

65 kms~i Mpc"^ 



3. Simulated data sets 



4. The x^-f't of cosmological parameters 



We use the numerical simulation package, SNOG, th e 
SuperNova Observation Calculator ( Goobar et al. , 2002a| ), 
to obtain simulated samples of the intrinsic dispersion and 
gravitational lensing effects of Type la SNe. The intrin- 
sic dispersion and measurement error is represented by a 
Gaussian distribution with am — 0.16 mag. Gravitational 
lensing effects are calculated by tracing the light between 
the source and the observer by sending it through a series 
of spherical cells in which the dark matter distribution can 
be specified. 

We model compact dark matter as point masses and 
smooth dark matter with the Navarro-Frenk- White fNFW; 



Navarro et al. (1997)] density profile using halo parame- 
ters, mass distributions and number densities as outlined 
in Bergstrom ct al. ( ^000 ). Note that the results obtained 
in this paper are not sensitive to the exact parameter- 
ization of the smooth halo profile. Also, the results are 
independent of the individual masses of the compact ob- 
jects as well as their clustering properties on galaxy scales 



( |Holz fc Wald| , |199^ ; [Bergstrom ct al.| , pOOOj ). Any eventual 
small subhalo structure in the dark matter halos does not 



act as a compact component ( [Mortsell et al.| , |2001a| ) 



Two different types of simulated data sets have been 
used for the analysis presented in this note. The first one, 
A, assumes that the number of observed Type la SNe 
is constant per comoving volume, which means that the 
number of events increases rapidly for higher redshifts. 
The B distribution instead assumes a constant number of 
SNe per redshift interval, i.e., a uniform z-distribution. 
Each data set consists of 2000 SNe in the redshift interval 
0.01 <z < 2.0. 

The number 2000 corresponds roughly to one 
year's data from the proposed satellite telescope 



Supernova/ Acceleration Probe [snap; Perlmutter ct al 
(200C)]. The exact redshift distribution of the Type la 



SNe to be followed by snap is the subject of several on- 
going science and instrumental optimization studies. The 
current anticipated distribution is roughly a constant rate 
per comoving volume for z < 1 and a uniform distribution 
at higher redshifts. In that respect is the A distribution to 
be regarded as a limiting case of maximal lensing effects 
due to the large number of high redshift SNe. 

To each of these data sets, 300 SNe between 0.04 < 
z < 0.08 have been added. These numbers correspond 



Making a x^-fit of {flM,^x,-M) using Eq. (|i|) (assuming 
w = —1) will for the data set B, describ ed in Sect. || 
result in the solid contour presented in Fig. |l(a)| if lensing 
effects are absent. If lensing dispersion is added, the fits of 
the cosmological parameters will be biased in accordance 
with the figure. All contours show the 1.51 a level, which 
corresponds to a '--^ 68 % confidence region of including the 
true value of the two parameters. 

The amount of bias depends on the matter distribu- 
tion responsible for the lensing effects. The dashed contour 
shows the result if only a new halo model is considered, 
while the dotted (dash-dotted) contour is based on a sim- 
ulation where the fraction of compact objects in the dark 
matter model is 20% (40%). 

Fig. l(b)| shows similar fits for the (f^M, w)-plane where 
a fiat universe is assumed {fix = 1 — ^m)- -M is still 
treated as a nuisance parameter, i.e., no prior knowledge 
is assumed. Also this figure is based on the data set B, 
but using the A distribution does not significantly alter 
the qualitative nature of the results. Gravitational lensing 
starts to become important for z > 1.0, and in the interval 
1.0 < z < 2.0, t he distributions A and B arc quite similiar. 

In Fig. 1(b) it is evident how important it is to consider 
lensing effects since a naive analysis of lensed data may in 
fact misleadingly exclude a cosmological constant as the 
explanation for the dark energy. 

The size of the bias, i.e., the systematic errors due to 
lensing in the above figures depends on the asymmetry of 
the magnitude distribution, which is described in Sect. ^ 
This in turn depends on the redshift for a given dark mat- 
ter model. The statistical error - the sizes of the ellipses - 
on the other hand only depend on the sample size so it is 
interesting to see where these two errors are comparable. 

In Fig. ^ the solid curve is the estimated statistical 
error of Qm plotted as a function of the sample size. It 
should be noted here that zero events means zero events 
from the high-z distribution, whereas the 300 events from 
the Supernova Factory simulations are always present. It 
is interesting to note in this figure that the full data set of 
2300 SNe is not necessary to get systematic errors compa- 
rable to the statistical error, but gravitational lensing has 
to be considered also for rather small data samples con- 
taining high redshift SNe (zj>l) if the fraction of compact 
objects is non-negligible. 
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(a) Fits in the (Hm, nx)-plane. 
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(b) Fits in tiie (JIm , ti')-plane assuming a fiat uni- 
verse. 

Fig. 1. Confidence contours showing the 1.51 ct level, i.e., 
the 68% level for two parameters, for three-parameter 
X^-fits where the third parameter Ai is treated as a nui- 
sance parameter. All fits are based on a simulated sample 
of 2300 SNe according to the data set B, but different 
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- ^ Bias for 40% CO. 
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Fig. 2. The evolution of the systematic error due to grav- 
itational lensing for different dark matter models and the 
statistical error (the solid curve) as a function of the sam- 
ple size. Note that zero SNe correspond to zero high- 
z events, i.e., the fits are only based on the 300 low-z 
events that are always added. The scatter in the lensing 
bias curves illustrates the stochastic nature of the effect: 
a handful of high-magnification supernovae may produce 
large effects in the fitted cosmological parameters. 



Universe, e.g., to determine the fraction of compact ob- 
jects like primordial black holes or MACHOs. 



As shown in Mortsell et al. (2001a), gravitational lens- 
ing effects are quite similar for different smooth dark mat- 
ter halo distributions, but very sensitive to the fraction of 
the matter density in compact objects, /p. Assuming that 
the major part of the matter density in the universe is 
in either smooth halos or compact objects, we can thus 
parametrize lensing effects with /p. 

Using SNOC, large data sets of synthetic SN obser- 
vations over a broad redshift range were created with a 
variable fraction of compact objects ranging from to 
40%, see Fig. ||. Using the Kolmogorov-Smirnov test, it 
was shown how the matter density in compact objects 
can be determined using SN data comparable to set A or 
B to < 5 % accuracy which we will assume to be the pre- 
dicted uncertainty in the paramete r /p, see also Metcali 



gravitc tional lensing models have been considered. The & Silk (1999) and 3eljak & Holz (1999). We now want 



dashed contours show the fit when a nfw halo model has 
been applied, and the dotted and dash-dotted contours 
represent different fractions of compact objects (co.) in 
the dark matter model. 



to parametrize the magnitude distributions obtained from 
the Monte-Carlo simulations to obtain approximately cor- 
rect pdf's to be used in the ML-analysis. 



5. Determining the lensing dispersion 

As we have seen, gravitational lensing may induce sizeable 
systematic effects when trying to determine cosmological 
parameters using the distance-redshift relation for stan- 
dard candle sources. A virtue of gravitational lensing is 
however that the distribution of luminosities can be used 
to obtain information on the matter distribution in the 



6. Parametrizing the pdf 

We parametrize the pdf's for different fractions of com- 
pact objects and redshifts with a Gaussian with a high 
magnification tail, i.e., 

/(to) — a-e 2^ , m > rric 

/(to) = a ■ e + b-\m\- lO' "", m < m^. (6) 
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Lensing+intrinaic dispersion [Amag] 



Fig. 3. Magnitude dispersion of data set B for % (full 
line), 20 % (dashed line) and 40 % (dotted line) compact 
objects using logarithmic scale. The bottom panel includes 
a Gaussian smearing, Um — 0.16 mag, due to intrinsic 
brightness differences between SNe and from measurement 
error. 

This distribution is motivated by the fact that the domi- 
nating intrinsic Gaussian dispersion will be shifted toward 
fainter values since a majority of lines-of-sight contain less 
matter than in a Friedmann-Lemaitre universe. The tail 
represents the effect of, e.g., galaxies lying close to some 
lines-of-sight causing high magnifications. 

We can set s — 1.5 and rric — for all cases and still 
have a reasonable y^. Since the fitted probability density 
functions are normalized to give an integrated probability 
of unity before performing the maximum likelihood anal- 
ysis, our parameters are toq, o and b/a. 

As an example, see Fig. ^, where Eq. (^ is fitted to a 
reference sample corresponding to 40 000 sources at z = 1 
in a universe with a 15 % mass fraction in compact objects. 
Here, mo = 0.01756, cr = 0.1653 and h/a = 0.02530 with 
xVndf = 86.44/68. 

For redshifts z < 0.5 (regardless of the value of /p) 
and 0.5 < z < 1 with low /p, one may ignore the high 
magnification tail (using a simple Gaussian) , and still have 
acceptable . Note however that the mean is still shifted 
from the zero value. For higher z and /p, we need the full 
pdf. E.g., for = 2 and 20 % compact objects, we get 
X^/ndf ~ 10 for a Gaussian and x^/ndf ~ 1 including the 
high magnification tail. 

7. Correcting for the bias 

There is apparently a great need to be able to handle 
gravitational lensing effects in the fits of the cosmologi- 
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Fig. 4. Fitting Eq. (|^) to the reference sample correspond- 
ing to 15 % compact objects. In this fit, a ~ 3597, mo = 
0.01756, cr = 0.1653 and h = 90.99. 

cal parameters. If the magnitude distribution presented 
in Sect. || is used in a ML-analysis the results presented in 
Fig. 1^ are obtained. The contours represent the ~ 68 % 
confidence region for one specific experimental realisa- 
tion.Q The solid contours are as previously the confidence 
region when no lensing effects have been considered while 
the dashed curve is the x^-?A of for a dark matter model of 
20% compact objects (cf. Fig. |l|). The dash-dotted curves 
show the ML-fits using the correct pdf, which is reducing 
the bias significantly as compared with the simple y^-fit. 
The dotted curves show the result of an ML-fit when frac- 
tions of 15 % and 25 % compact objects are assumed when 
the true value is 20%. These contours represent the ex- 
pected systematic bias when under- or overestimating the 
fraction of dark matter in compact objects with 5 %. The 
results are a big improvement over the x^-analysis. Thus, 
it should be possible to reduce the bias due to lensing in 
a fairly effective way, even with the simple approach used 
in this note. 

8. Summary and discussion 

Neglecting gravitational lensing effects may cause system- 
atic errors when determining cosmological parameters us- 
ing the distance-redshift relation for Type la SNe. E.g., if 
the universe contains a large fraction of compact objects, 
the matter density SIm can be severely underestimated or 
a cosmological constant may be wrongly excluded if one 
assumes a Gaussian distribution of magnitudes. Using an 
approximately correct pdf for the SNe magnitudes, this 
bias can be reduced significantly. 

In this note, we show that lensing effects can be 
parametrized in a simple way by the fraction of compact 
objects, /p, and the redshift. In general, lensing effects 
are bigger for larger values of /p and higher redshifts. If 

^ A different realisation would give contours of the same size 
but these may not have the same centres. They should however 
also include the true value which gives a limit on the how much 
of a difference there could be between different realisations. 
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/p can be determined with < 5 % accuracy, the bias in the 
cosmological parameter determination can be effectively 
reduced as compared to the case where the magnitudes 
are assumed to be Gaussian distributed. For future sur- 
veys, e.g., by the proposed SNAP satellite aiming at doing 
precision cosmology with Type la SNe, this reduction is 
essential to get systematic errors comparable to the sta- 
tistical errors. 

For the simple analysis performed in this note, we 
have first determined /p to be able to choose the cor- 
rect pdf for the subsequent parameter determination. 
Another possibility is to make a simultaneous fit of, e.g., 
(/p, r^M, ^x,.M) which will be the aim of some forthcom- 
ing work. 
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(a) The fits in the (SIm, f2js:)-plane. 
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(b) The fits in the (JIm, '!i')-plane, assuming a flat 
universe. 



Fig. 5. Both plots show the result of three-parameter fits 
where A4 has been fitted as a nuisance parameter. The 
solid contours are the x^"fit result when lensing effects 
are absent and the dashed contours are the x^-fit for a 
dark matter model with 20% compact objects and 80% 
parametrized as nfw if lensing effects are neglected in 
the fit (cf. Fig. [|). The dash-dotted contours show the 
fit result on the same sample when the parametrized pdf 
presented in Sect. ^ is used and the correct value of 20% 
compact objects are used. The dotted contours show the 
result when an incorrect value (±5 %) of the fraction of 
compact objects is assumed. Note that the confidence re- 
gions represent the outcome of one specific experimental 
realisation. 



